Chapter 1 

Carbon nanotubes: 
Nonlinear high-Q 
resonators with strong 
coupling to single-electron 
tunneling 

Harold B. Meerwaldt, Gary A. Steele, and Herre S.J. van der Zant 

1.1 Introduction 

A carbon nanotube (CNT) is a remarkable material, which consists of only 
carbon atoms and can be thought of as a graphene sheet rolled up into a cylinder 
with the ends capped off with a buckyball sliced in half. This is called a single- 
walled CNT and, depending on how the graphene sheet is rolled up, CNTs are 
either semiconducting or metallic. At low temperatures, CNTs contacted by 
two electrodes become quantum dots, and Coulomb blockade and single-electron 
tunneling occurs. When more graphene sheets are wrapped up concentrically a 
multiwalled CNT is formed. 

CNTs are fabricated using various methods. These include arc discharge 
PP, laser ablation [2], chemical vapor deposition (CVD) |3], and the use of high 
pressure carbon monoxide [4]. There are two ways to make doubly-clamped 
CNT resonators. The first method starts with growing CNTs on an oxidized 
silicon wafer and then undcrctching them with a trench so that they become 
suspended. This method exposes the CNT to resist, the electron beam, and 
acid, and may give rise to defects in and residues on the CNT. In the ultraclean 
method [5] , processing on the CNT is avoided by first making the contacts and 
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the trench and in the final step growing the CNTs. A suspended CNT res- 
onator can be fabricated without defects, thus reducing damping which results 
in quality factors for the flexural modes above 100,000 [6]. 

Having a radius that can be as small as a nanometer and lengths of several 
micrometers, CNTs have aspect ratios in the thousands, on a scale and purity 
that is difficult to achieve in resonators fabricated top-down from silicon, SiN, 
or other materials. The enormous aspect ratio makes it easy to excite the 
CNT such that the displacement of the flexural mode is similar in magnitude 
to the radius; the CNT thus behaves as a thin, narrow beam resonator. The 
combination of low damping and high aspect ratio paves the way for nonlinear 
effects. 

There arc two other remarkable properties of CNTs: Their Young's mod- 
ulus of 1.3 TPa makes them an incredibly stiff material, and combining this 
with a density of only 1350 kg/m 3 , CNT resonators reach frequencies of sev- 
eral hundreds of megahertz and even gigahertz. At millikelvin temperatures, 
with Hujq 3> kt,T, thermal phonons can no longer excite the flexural motion. 
This makes CNTs perfect candidates to observe quantum (nonlinear) effects in 
mechanical resonators [7]. 

This chapter is built up as follows. In the next section various ways of detect- 
ing the motion of CNT resonators are discussed. In the third section we show 
how single-electron tunneling gives rise to frequency softening and damping. 
Section four contains Duffing-like nonlincaritics caused by geometry, electro- 
statics, and single-electron tunneling. In section five we look at nonlincaritics 
from parametric driving and mode coupling to other dimensions of oscillation 
and to higher modes. Table [TTT1 gives an overview of the symbols used in this 
chapter for variables and constants along with typical values and example device 
parameters. 

1.2 Detecting the motion of a carbon nanotube 
resonator 

In this section, we discuss different ways to detect the flexural motion of CNT 
resonators. Other vibrational modes such as breathing and stretching modes 
are also present in CNT resonators, but non-flexural modes are difficult to ac- 
tuate and detect. They are therefore not discussed in this chapter. Because a 
CNT is easily perturbed, detection has a large influence on its motion. When 
discussing the nonlinear dynamics of CNT resonators it is therefore important 
to have an overview of the different detection methods. Optical detection, which 
is used frequently in top-down devices, is not an option because of the small 
cross-section of a CNT. Until now, detection has therefore been done through 
microscopy or electronically. We summarize five different methods to observe its 
motion: scanning force microscopy, transmission electron microscopy, field emis- 
sion microscopy, the mixing technique, which involves frequency down-mixing 
the flexural motion from several hundreds of megahertz to several kilohertz, 
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General properties of single-walled CNTs 



E 


Young's modulus 


1.3 TPa 


P 


mass density 


1350 kg/m 3 


Example CNT 


L 


length 


800 nm 


r 


radius 


1.5 nm 


hg 


distance to gate 


230 nm 




capacitance to gate 


10 aF 


Ctot 


total capacitance 


20 aF 


Q 


quality factor 


100,000 


loq/2t: 


resonance frequency 


140 MHz 


I 


second moment of inertia 


4.0 nm 4 


A 


cross-sectional area 


7.1 nm 2 


Variables 


u(z, t) 


displacement 




A{t) 


slowly-changing amplitude 




i{z) 


mode shape 





Table 1.1: General properties of a single- walled CNT, device parameters for a 
CNT used as an example in the subsequent sections, and variables. 

and the rectification technique, where the amplitude of the flexural motion is 
obtained by measuring the dc current flowing through the CNT. The last two 
(electronic) methods are self-detecting because the CNT is both the object stud- 
ied and the detector. 

Using scanning force microscopy (SFM) jS], it is possible to spatially image 
the different mode shapes in a CNT resonator. Figure 11.2.11 shows an SFM in 
tapping mode with the tip positioned above a doubly clamped CNT. The CNT 
is actuated by a nearby gate at a frequency uj. Since the bandwidth of the SFM 
is not high enough to measure the oscillations of the CNT at the resonance 
frequency, an amplitude-modulated actuation signal is used for the voltage on 
the gate with a frequency of ui m od- The tip cannot follow the high-frequency 
CNT vibrations, but it can follow the envelope of the oscillation. On resonance, 
the amplitude-modulated actuation results in a fast-oscillating displacement 
envelope u with frequency uimod (see figure ll.2.1[) . The modulation frequency 
can be conveniently chosen such that the signal can be measured with a lock-in 
amplifier. For higher sensitivity, u> mo d is matched to the first eigenmode of the 
SFM tip. The amplitude of oscillation at different locations on the CNT provides 
a time averaged image of the shape of the driven mode. With transmission 
electron microscopy (TEM) [9], a fiber of CNT cantilevers is connected to a 
gold wire, which is placed between the electron gun and the viewing system, as 
shown in figure fl~2~2] By applying a voltage to the CNTs, the ones that are not 
perpendicular to the grounded counter electrode are attracted to it. To drive 
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Figure 1.2.1: Scanning force microscopy. A suspended doubly- clamped CNT is 
brought into motion by a nearby gate electrode (left). The amplitude-modulated 
motion (right) at different parts of the CNT is read out at 0J mo d by a tip in 
tapping mode. 

the CNTs, an alternating voltage is applied to them and the resulting envelope 
of the motion is observed in the TEM image. As with the SFM method, the 
low bandwidth of TEM imaging allows visualizing only of the mode shape of 
the vibrations. 

Reports on these two methods have not given indications of nonlinearity. In 
the experiment using SFM quality factors in the order of 10 and the disturbance 
due to the tapping tip may have prevented nonlinearity from being observed. In 
the experiment with TEM, the high aspect ratio of the CNT makes it difficult 
for the CNT cantilever to bend to such an extent that geometric nonlinearity 
becomes relevant [10]. Even though the amplitude of oscillation is far larger 
than the CNT radius, tension induced by bending is much smaller in CNT 
cantilevers than in doubly-clamped CNTs. 

With field-emission microscopy (FEM) , temporal and spatial information on 
the flexural motion is obtained. Figure 11.2.31 shows a CNT cantilever to which 
a few hundred volts is applied. A field-emission tunnel current now flows to the 
anode which is placed up to a few hundred nm away. In one case jll) , the CNT 
was excited by radio waves from a nearby antenna. At the mechanical resonance 
frequency the distance between the tip of the CNT, which acts as a cathode, and 
the anode oscillates. This causes the tunnel resistance and therefore the current 
flowing through the circuit to oscillate. In this way, the CNT acts as a rectifier; 
the current is highest when the CNT is in its equilibrium position and becomes 
lower, but never negative, when the CNT moves away from this position in 
either direction. As a consequence, the time-averaged current decreases when 
the CNT is set in motion. Thus, the ac motion of the the CNT is probed by 
a dc measurement of the averaged current. In another case the CNT was 
positioned near a phosphor screen so that the averaged 2-dimensional image of 
emitted electrons was recorded by a 25 Hz video camera. 

It should be noted that the three techniques discussed up to now operate at 
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electron beam 




Figure 1.2.2: Transmission electron microscopy. A suspended CNT cantilever 
is driven by applying a voltage Vcnt to the CNT with respect to a grounded 
counter electrode. Electrons from an incident electron beam are not transmitted 
through the CNT and as result an image of the flexural motion is obtained. 

V RF cos(wt) — f 



Figure 1.2.3: Field emission microscopy. A suspended CNT cantilever is excited 
by radio waves from a nearby antenna. Current tunnels from the CNT to the 
anode and oscillates with the motion of the CNT resonator. 

room temperature. At low temperatures, mechanical properties such as damping 
change dramatically. We now discuss two electronic methods that are used at 
low temperature, which we refer to as the mixing technique and the rectification 
technique. 

In measuring the high-frequency response of mechanical resonators, one of- 
ten employs electrical frequency-mixing techniques. Here, a nonlinear element 
mixes two different high-frequency signals to their difference frequency, which is 
considerably lower and is measured more easily. In the case of a CNT mechani- 
cal resonator, the CNT itself is used as the nonlinear element. To measure CNTs 
using the mixing technique |13[ I14[ I15[ 116) , the working principle is modeled as 
follows. A three terminal setup is used, as shown in figure 11.2.41 The current 
/ = GV flowing through the CNT, where G is the conductance of the CNT 
and V is the bias voltage over the CNT, depends on the gate induced charge 
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<lc = CgVg, where C g is the capacitance to the gate and V g is the gate voltage. 
When the CNT moves away from its equilibrium position, the capacitance to the 
gate changes, giving q™ ech = C g c (t)V g . Also, when an oscillating signal V g c (t) 
is applied to the gate, the gate induced charge changes as: q^ irect — C g V g c (t). 
Putting these two contributions together gives the time-dependent conductance: 

^)=§X v °°' {t)+ W v ')- 

The two contributions both oscillate at a frequency U) but can in general have 
a different phase. Note that the frequency of the mechanical oscillation of the 
CNT equals the frequency of the driving voltage applied to the gate electrode. 
When also the bias voltage is alternated, but now at a frequency u + Aw, 
both the mechanical oscillations and the direct oscillations are mixed given 
that / = G ac (t)V ac (t). One mixing component occurs at Aw, which can be 
conveniently chosen at a few kHz. This signal is then detected using a lock- in 
amplifier. At high driving, nonlinear effects become visible [131 115] . 

Measurement of the flexural motion of a CNT resonator has also been per- 
formed using the mixing technique with a frequency-modulated voltage applied 
to the source electrode [17]. The frequency modulation of the source voltage 
causes sidebands in the spectrum around the drive frequency. They both mix 
with the mechanical resonance and result in a signal at the modulation fre- 
quency. In this case no signal needs to be applied to the gate as the CNT is 
actuated by setting the carrier frequency of the source voltage to the mechanical 
resonance frequency. 
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Figure 1.2.4: Detecting the CNT's motion using the mixing technique. A CNT 
is doubly clamped between a source and a drain electrode suspended above 
and driven by a gate electrode. The flexural motion of the CNT mixes the 
bias voltage at w + Au and the gate voltage at u to a current at Au. When 
frequency-modulation of the source voltage is performed, no ac voltage needs to 
be applied to the gate electrode. 

Another technique that resembles the mixing technique transduces the me- 
chanical oscillations of the CNT resonator at several hundred megahertz into a 
dc current [6], in this way using the CNT resonator as a rectifier. The advan- 
tage of this technique is that the amplitude of the mechanical motion is given 
by a change in the dc current flowing through the CNT. As with the mixing 
technique, the use of the rectification technique avoids the difficulty of getting 
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small, high-frequency signals out of a setup at millikelvin temperatures and 
minimizing the cross-talk from the actuation onto the measurement signal. 

With the rectification technique, the working principle is as follows. The 
CNT is suspended between source and drain electrodes above a gate electrode. 
At low temperatures, the CNT acts as a suspended quantum dot in which charg- 
ing effects dominate transport. It is actuated using a nearby antenna, which 
sends out an oscillating electric field. When the CNT oscillates the distance to 
the gate changes, and therefore the capacitance Cg C (t) = -^f-u ac cos(ivt) , where 
u is the displacement of the CNT. Since the only relevant quantity is the charge 
induced by the gate, it is equivalent to say that the gate voltage changes as 

C ac (t) 

Vg C (t) = ~ c — ^9' Expanding the current through the CNT quantum dot with 
respect to gate voltage one obtains: 

Id 2 / fV g dC g \ 2 (\ 1 

Us, 



- + -cos(2^) . (1.2.2) 



2dV g 2 \C g du ) ac \2 2 

The high resistance of CNTs combined with a high capacitance to ground often 
means that the bandwidth is too low to measure the oscillations at u> and 2uj 
directly. Omitting these oscillating terms gives for the change in the dc current 
AI due to mechanical motion: 

Important to note is that the change in current is averaged over time and that 
it is proportional to the amplitude squared. 

Equation ll.2.31 shows that, in order to observe a large mechanical signal, it 
is required that d 2 I/dV 2 be high, as is the case in the Coulomb blockade regime 
[T8l [T9] . When measured on a Coulomb peak the second derivative is negative 
and therefore the change in current as well. In the Coulomb valley both the 
second derivative and the change in current are positive. Typically, d 2 I/dV g 2 is 
several iiA/V 2 . On the inflection point of the Coulomb peak d 2 I/dV 2 = and 
there is no change in current. This is in contrast to the mixing technique which, 
as shown in equation 11.2.11 has the highest signal on the inflection point and 
the lowest on the Coulomb peak or in the valleys. 



1.3 Linear dynamics of a carbon nanotube res- 
onator strongly coupled to a quantum dot 

Because of their high quality factor and therefore narrow linewidths, ultraclean 
CNT resonators allow for the detection of small changes in resonance frequency. 
In this section we describe the linear response of CNT resonators measured with 
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Figure 1.2.5: Detecting the CNT's motion using the rectification technique. 
A CNT is doubly clamped between a source and a drain electrode suspended 
above a gate electrode. The CNT is driven by a high-frequency signal applied 
to an antenna at approximately 1 cm from the CNT. The single-electron tunnel 
current through the CNT quantum dot contains information on the oscillation 
of the CNT at dc due to a rectification effect from Coulomb blockade. 



the rectification technique. These measurements are performed at low temper- 
atures; the CNT then acts as a quantum dot [20l [211 E2], and single-electron 
tunneling is the dominant transport mechanism. We consider the backaction 
force caused by tunneling electrons and how this stochastic force gives rise to 
dips in the resonance frequency and increased damping. 

Figure 11.3.11 shows Coulomb oscillations (top) , and the change in resonance 
frequency with gate voltage (bottom) at negative gate voltages. Three features 
are worth mentioning: i) in Coulomb blockade, there is a fixed charge on the 
CNT and the resonance frequency-gate voltage curve approaches a linear de- 
pendence, denoted by the dashed lines; ii) from one charge state to the next, 
there is an offset between the sloped lines; and iii) at the position of Coulomb 
peaks corresponding to a transition between charge states, there are resonance 
frequency dips arising from a softening-spring behaviour. 

The slope in resonance frequency is explained as follows: in Coulomb block- 
ade the charge on the CNT quantum dot is constant. If a more negative gate 
voltage is applied, the force acting on the CNT increases, which pulls the CNT 
more towards the gate. The tension caused by this pulling leads to the first 
feature: the resonance frequency increases with gate voltage. 

As the gate voltage becomes more negative, more and more electrons are 
extracted from the CNT quantum dot, leaving behind holes. From one charge 
state to the next, there is a difference in charge of one electron. The tension 
arising from an extra electron causes the resonance frequency to change in a 
discrete manner, leading to the second feature: an offset in resonance frequency 
between two charge states. Figure [T.3.11 shows a frequency offset of 0.5 MHz, 
which is more than 100 times the lincwidth. 

To explain the frequency dips arising from single-electron tunneling, we look 
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Figure 1.3.1: Top: Coulomb oscillations of the current through a CNT quantum 
dot; in Coulomb blockade current is zero and the indicated (negative) number of 
electrons on the CNT is fixed. Between charge states a single-electron tunneling 
current is visible as Coulomb peaks. Bottom: Gate dependence of resonance 
frequency (from [13]); the features are indicated by the encircled numbers: (1) 
a slope caused by tension at fixed charge, (2) a frequency offset caused by 
a tension difference between charge states, and (3) frequency dips caused by 
single-electron tunneling. 



at the electrostatic force acting on the CNT: 

where dC g /du is the derivative of the capacitance to the gate with respect to 
displacement and V g is the gate voltage. The voltage on the CNT, Vcnt, is 
given by: 

V CNT = C ' V « + qcNT , (1.3.2) 

CcNT 

where qcNT is the charge on the CNT and Ccnt = C g + C s + Cd is the total 
capacitance of the CNT to gate, source, and drain. The charge on the CNT is 
given by: 

qcNT = -eN(q c ), (1.3.3) 

where e is the elementary charge, and N(q c ) is the number of electrons on 
the CNT, which is negative in figure 11.3.11 The charge and voltage on the 
CNT are shown in figure 11.3.21 The gate induced charge q c is the charge that 
would be on the CNT if there were no Coulomb blockade, and in the case of 
antisymmetrically applied voltages (V s = — Vd) and equal capacitances (C s = 
Cd) for source (s) and drain (d), it is given by: 

q C = CgVg. (1.3.4) 

As the CNT undergoes a charge transition, electrons can hop on and off the 
CNT quantum dot. Whereas in Coulomb blockade there was no possibility 
to screen the voltage on the gate, in a charge transition the CNT overscreens 
the gate voltage, so that on average the voltage on the CNT remains zero. It is 
this overscreening which causes a reduction in tension and gives rise to the third 
feature: dips in the resonance frequency due to single-electron tunneling. Figure 
11.3.11 shows that single-electron tunneling can lead to a dip in the resonance 
frequency by as much as 2 MHz. The change in resonant frequency due to 
single-electron tunnelling over this small gate range is an order of magnitude 
larger than that from the gate induced mechanical tension, which is described 
in section [1. 4. 21 

By considering the electrostatic force, the change in spring constant and 
thereby the change in resonance frequency due to single-electron tunneling can 
be calculated: 

a; dF C NT /t r t r ^dCgdVcNT 
kksET = j ~ (V g - Vdot) , , 

du du du 

w VM-Vcst) f^Vfl-e^Ws) 
Ccnt \ du J \ dq c J 

where we neglect slowly varying terms with d 2 C g /du 2 . 

The average number of electrons, or average occupancy (JV), can be non- 
integer (dashed line in figure [TT3T2]) . At fixed — qcNT> d(N)/dq c equals zero and 
the factor V g (V g — Vcnt) represents the tension causing the first two features, 
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Figure 1.3.2: Left: The charge on the CNT increases with integer steps as the 
gate induced charge is increased (solid line). At finite temperature or finite 
tunnel coupling the charge transitions are smoothed out (dashed line). Right: 
The voltage on the CNT increases with increased gate induced charge as Vc nt = 
CgVg/CcNT until an extra charge enters the CNT, at which point the CNT 
voltage decreases with e/CcNT (solid line). At finite temperature or finite 
tunnel coupling this voltage jump is smoothed again (dashed line). 

the frequency slope and the frequency offset. At a charge transition, d(N)/dq c is 
positive, and this gives rise to softening-spring behaviour in equation 1 1.3. 51 The 
sharper the charge transition is, the deeper the frequency dip is. This follows 
from the fact that when d(N)/dq c is higher, the average number of charges on 
the CNT changes more as the CNT oscillates towards and away from the gate, 
in this way giving rise to a larger change in force with change in displacement, 
dFcNT I 'du, and therefore a larger change in spring constant. 

Using A A; = 2mu>oAujQ we arrive at the change in resonance frequency due 
to single-electron tunneling: 



In addition to a change in resonance frequency, electrons hopping on and off the 
CNT quantum dot can also cause damping of the mechanical motion. We look 
at the limit when there are many tunneling events per mechanical oscillation, 
T to t ^ ^o- As derived in [21], the total damping has an intrinsic contribu- 
tion and a contribution due to the stochastic backaction force associated with 
tunneling electrons: 



Vg(Vg-VcNT) ( dC g 

2mujoCcNT \ du 




2 




(1.3.6) 




(1.3.7) 
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where Qtot is the total quality factor, Qmt is the quality factor in the absence 
of tunneling electrons, is the rate of electrons tunneling onto (off) the CNT, 
and r tot = r + + T~ is the total tunneling rate. The stochastic force the CNT 
experiences, F s toch = F(N + 1) — F(N), is the difference between the force expe- 
rienced at N and N+l electrons. Assuming that both source and drain voltage 
are much smaller than the gate voltage, V s , Va -C Vg, and taking only into 
account the dc force experienced by the electrons as the CNT oscillates (which 
is valid as long as V g ac <C V g dc ), reference [24] gives the following expression for 
the stochastic force: 

F s toch = 7^-^ (2e(C L + C R )V B + e 2 (2N + 1)) . (1.3.8) 

In this limit, the stochastic force on the CNT depends on the distance of the 
CNT to the gate and consequently on the average number of charges on the 
CNT, both of which change as the CNT oscillates. To obtain F stoc h for holes, 
the substitution e — > —e must be made in equation ll .3.81 with N now the number 
of holes. Using 

1 /v r 

(N) = — — '— - — = — v -, (1.3.9) 

v ' i/r+ + i/r_ Ttot 1 ; 

this leads to the following expression for the quality factor: 

Q**=( 7 ±- + /V 2 _ ffiV (MCl + C R )Vt + e\2N + 1)) [1.3.10) 

\Qint mu r to tCf ot C g \ du ) V y > dV g J 

Figure 11.3.31 illustrates how a charge transition from iV to N + 1 gives rise to a 
decrease in the quality factor. Similar to the frequency dips discussed earlier, 
damping is highest when d(N) /dV g is highest. At this point, the mechanical 
oscillation causes the largest change in the average occupancy and consequently 
the influence of the backaction is the most pronounced. Physically the damping 
is caused by energy dependent tunneling. When the CNT is oscillating, the 
change in capacitance to the gate makes the chemical potential of the CNT 
oscillate. Even though many electrons tunnel onto and off the CNT during 
one mechanical oscillation, on average more tunnel on when the energy level is 
lower, and more tunnel off when the energy level is higher. This asymmetry is 
imposed by the Fermi energy of the source and drain. Figure H .3.41 shows that, 
when the chemical potential of the CNT is lower than the Fermi level of the 
source or drain, electrons are allowed to tunnel onto the CNT but not off the 
CNT. Conversely, when the chemical potential of the CNT is higher than the 
Fermi level of the source or drain, electrons are allowed to tunnel off but there 
are no electrons to tunnel onto the CNT. Because of this asymmetry, electrons 
will tunnel on when the level is pulled below the Fermi energy by the motion, 
and tunnel out when the level is pushed above the Fermi energy by the motion, 
thereby extracting energy from the mechanical mode, damping it's motion. Note 
that this energy is carried away by hot electrons above the Fermi energy in the 
leads, where it is later dissipated when these electrons relax. We finally note 
that the gate-dependent damping, as illustrated by the black curve in figure 
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Figure 1.3.3: Noninteger part of the average occupancy (left axis) at a charge 
transition from TV to TV + 1 electrons and the corresponding change in quality 
factor (right axis) through single electron tunneling for the example CNT with 
20 electrons at V g = 3V assuming r tot = 10 13 Hz. The maximum of the Coulomb 
peak is located at mV of gate voltage offset. 

11.3.31 has been observed in experiments. The change in Q can be dramatic: 
from 10 5 in the Coulomb valleys to several hundreds near the Coulomb peak. 
At present we are performing a more accurate comparison with the theoretical 
calculations [25] . 

1.4 Carbon nanotube resonators in the nonlin- 
ear regime 

Hooke's Law states that, when an elastic object is brought out of its equilibrium 
position, the restoring force is linearly dependent on the displacement. As the 
applied force and resulting displacement increase, Hooke's Law is no longer valid 
and nonlinear effects become important. The first correction to Hooke's Law 
is made by adding a restoring force term that is quadratic in displacement. 
The quadratic term leads to a renormalization of the cubic term [26j . In many 
cases, this renormalization is a small effect, and hence we have neglected it in 
this analysis. In this section we focus on the next term of the restoring force, 
which is cubic in displacement. The cubic term is called the Duffing term, 
and the resulting equation of motion the Duffing equation (see also [26]). The 
restoring force is now given by F r = — k^u — mau 3 , were a is the nonlinearity 
or Duffing parameter with units of N/(kg m 3 ). The nonlinearity can also be 
seen as a correction to the spring constant. Writing the restoring force as 
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Figure 1.3.4: Energy diagram (exaggerated): Due to the asymmetry imposed 
by the Fermi level of the source and drain, electrons are allowed to tunnel onto 
the CNT but not off the CNT when the chemical potential of the CNT is lower 
than the Fermi level of the leads. When the chemical potential of the CNT is 
higher than the Fermi level of the leads, electrons are allowed to tunnel off the 
CNT but there are no electrons to tunnel onto the CNT. 



F r = —ku = — (fen + motu )u, the nonlinearity parameter is given by: 

2m du 2 

When a is positive, the Duffing term causes, as drive force is increased, 
the resonance peak to tilt towards higher frequencies, resulting in a stiffening- 
spring behaviour. Conversely, when a is negative and the resonance peak tilts 
towards lower frequencies, it results in a softening-spring behaviour. As the 
applied force is increased further, the tilting of the peak causes bistability: for 
one drive frequency there are two metastable states of oscillation, one with a 
low amplitude and the other with a high amplitude of oscillation [7] . 

In this section we first describe how the measured lineshape changes as the 
drive power is increased and at what power the lineshape becomes nonlinear. 
Next, we describe three causes for nonlinearity in CNT mechanical resonators: 
geometric nonlinearity caused by mechanical tension, electrostatic nonlinearity, 
and nonlinearity due to single-electron tunneling. We calculate the correspond- 
ing a values for the example CNT described in table 11.11 



1.4.1 Nonlinear response of carbon nanotube resonators 

Figure 11.4.11 shows a measurement of the lineshape of a 600 nm long CNT res- 
onator detected by using the rectification technique (see section [O]). At a low 
power of 
-43 dBm, the CNT is driven to oscillate at a small amplitude and the noise 
is too high to observe the motion. As the drive power is increased to -40 dBm, a 
Lorentzian lineshape appears in the current when the drive frequency is near the 
mechanical resonance frequency of the CNT. At this drive power, the restoring 
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force is given by Hooke's Law and can be approximated to be linear with dis- 
placement. When the drive power is increased further to -38 dBm, the restoring 
force is nonlinear with respect to displacement but the amplitude is still single- 
valued, resulting in a sharkfin-like appearance of the lineshape. At a higher 
power of -35 dBm in figure 11.4.11 bistability and hysteresis occur. Bistability 
appears as for certain frequencies there is both a low- and a high-amplitude state 
in which the resonator can oscillate. The amplitude of the motion depends on 
the direction in which the drive frequency is changed, displaying a hysteretic be- 
haviour when the drive frequency is swept from high to low, or from low to high. 
The bottom two panels of figure ll~4.il show nonlincarity in the form of a Duffing 
sharkfin-like lineshape. The sharkfin lineshape points towards higher frequency, 
indicating a stiffening spring behaviour. In many applications "230 ["TS] , the range 
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Figure 1.4.1: Frequency responses of the rectified current (see previous section) 
showing the mechanical signal drowned in noise (-43 dBm), a linear response 
with a Lorentzian lineshape (-40 dBm), the onset of nonlincarity (-38 dBm), 
and nonlinearity with bistability and hysteresis (-35 dBm), of a CNT across a 
600 nm wide trench, 285 nm above a gate, measured at 20 mK. 



of powers over which the mechanical resonance displays a detectable linear re- 
sponse is an important parameter. This range of powers is often referred to 
as the "dynamic range" of the mechanical oscillator. The smallest useful drive 
power is the one for which the measured signal matches the measured noise, 
and the largest useful drive power is the one at the onset of nonlinearity. 

To estimate an upper bound for the dynamic range, we follow the analysis 
of |29| and only consider the thermomechanical contribution to the noise, leav- 
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ing out other intrinsic and extrinsic sources of noise. The spectral density of 
displacement noise at the resonance frequency is: 

S, - (1.4.2) 

where kb is Boltzmann's constant, T is temperature of the resonator, Q is the 
quality factor, m is the total mass of the resonator, and ujq is the resonance 
frequency For the example CNT in table 11.11 this gives a displacement noise of 
y/SZ = 0.15 pm/VHz. The root-mean-square amplitude of oscillation at which 
nonlinearity sets in, A c , is the amplitude at the lowest drive power at which an 
infinitesimal change in frequency results in an infinite change in amplitude. The 
analysis of nonlinear response curves and the onset of nonlinearity can be found 
in [30]. According to [31] the critical amplitude, which is not at the resonance 
frequency, is given by: 

. 2 ■ 0.83 /2V3 7 
A= = ^3- W0 V^Q' (1.4.3) 

where a is the nonlinearity parameter which will be discussed in the next section. 

The dynamic range is defined as the ratio of the 1 dB compression point to 
the noise amplitude at the resonance frequency: 



DR(dB) = 20 log I 



' 0.745 Ac 
0.745- 0.83 qjg / Z^mup \ 



20 3Q < L «> 

where A/ is the measurement bandwidth. The 1 dB compression point occurs 
when the peak amplitude is 1 dB lower than what would be expected for a 
purely linearly response and has a value of 0.745„4 C . Equation 11.4.41 shows that 
the low mass and the high nonlinearity parameter for CNT resonators give rise 
to a small dynamic range compared to top-down fabricated nanomcchanical 
resonators. 



1.4.2 Geometric nonlinearity 

In this subsection, it is described how geometric nonlinearities in the motion of 
CNT resonators arise. When the amplitude of oscillation becomes larger, the 
CNT experiences significantly more tension at the extremal points of oscillation 
than at the equilibrium position. This tension thus depends on displacement 
and is one contribution to the geometric nonlinearity. Furthermore, the force 
induced by the gate voltage causes tension in the CNT, which again depends 
on displacement, being another contribution to the geometric nonlinearity. 

To derive the expressions describing geometric nonlinearity we start with 
the theoretical framework for a driven damped beam [32j . The Euler-Bernouilli 
beam theory describes the motion of continuous bodies subjected to external 
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forces. Even though the diameter of a CNT consists of only around 10 atoms, 
continuum mechanics is nonetheless a good description of the mechanical prop- 
erties of a CNT resonator [33]. We take the coordinate system as shown in 
figure [T74.2I The CNT axis is positioned in the z-direction and the oscillation 
takes (at least for now) place in the x-direction. The displacement u of the CNT 



x CNT 



source 






drain 




rh 
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insulator 
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Figure 1.4.2: Coordinate axes: the CNT resonator of length L is suspended 
between the source and drain electrodes at a distance of h g above the gate. 
The CNT is aligned in the z-direction. The displacement u takes place in the 
x-direction. 



resonator as a function of time is governed by the following equation of motion: 

, d 2 u du ^d 4 u m <9 2 u 
pA W + ^ +D d^- T d^ = f - 

The first term describes inertia, with p being the density and A the area of 
the cross section. The second term describes damping with 77 = pAujo/Q the 
damping rate per unit of length. The third term is the restoring force due 
to the bending rigidity D = EI /(I - v 2 ), with E = 1.3 TPa [34] and v = 
0.2 [35] the Young's modulus and Poisson's ratio respectively, and the second 
(bending) moment of inertia /. The CNT's Poisson's ratio is neglected and 
the denominator is set to unity. We approximate the CNT resonator as a solid 
cylinder with a second moment of inertia I = 7rr 4 /4. The fourth term describes 
the restoring force due to the tension T and / describes an external force per 
unit length. 

There are four boundary conditions determining the modeshapes u(z). The 
CNT is doubly clamped, leading to it(0) = and u(L) = 0. Furthermore, the 
bending rigidity of the CNT prescribes the slope of the CNT to be continuous 
at the clamping points, as the CNT goes from lying on the contact to being 
suspended, resulting in du/dz(0) = and du/dz{L) = 0. 

To solve equation 11.4.51 we first separate all contributions into a static/dc 
and a dynamic/ac component. For the position dependent displacement of the 
CNT, u(z), and for the gate voltage, V g (z), this gives: 

u(z) = Udc(z) + Uac(z) cos(wt), (1.4.6) 
V g (z) = V g dc (z) + V g ac (z)cos(ujt). (1.4.7) 
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The tension felt by the CNT is then given by: 



EA f L fdu^ 2 



z 




rp . EA p fdu dc y EA I f du dc du a 

= T dc + T ac cos(ut), (1.4.8) 

with To the residual tension present at zero gate voltage. The uniform force per 
unit length on the CNT due to the gate voltage is: 

= fdc(z) + f ac (z)cOs(uit). 

Here dc g /du is the derivative of the capacitance per unit length with respect 
to displacement. We arrive at two equations; one for the static and one for the 
dynamic displacement: 

D ^- TdC ^f =fdC > (1A10) 

2 a , n d 4 tt ac d 2 u ac d 2 u dc 

-UJ pAu ac + lUJ?]U ac + D—— T dc — — T ac g 9 = Jac- (1.4.11) 

az 4 oz z oz z 

It is useful to introduce scaled variables denoted with primes: 



z 



rp/ T dc ,, 



L ' dc,ac D 



V^-T. f=*fi. = & (1-4-12) 

The equation for the dynamic displacement now becomes: 



u u 



, d 2 u' ac d^_ d 



7? U ac + U— — l dc l ac - L ac . \IA.16) 



The ac displacement is separated into an amplitude A(t) changing slowly with 
respect to the oscillation, and a spatial part £,(z): 

u ac (z,t) = A n (t)£ n (z), (1.4.14) 

where the subscript n denotes the mode of oscillation. The cigenfunctions of 
the spatial part of equation II .4. 131 give the modeshapes: 

CUZ ' ] = { D &± ~ T ' d ^) 61 ~ T ^n}^j£ = "n{T' dc , u' dc ) ■ Uz\ 

(1.4.15) 
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where the eigenfrequencies are given by: 



w' n =ftJl + '%?, (1-4.16) 

with p n a numerical factor. The first four /?„ are found to be 4.73004, 7.8532, 
10.9956, and 14.1372. Like a guitar string, the resonance frequency of the CNT 
increases with increased tension. Also, as tension is increased the resonance 
frequencies converge towards a harmonic spectrum. 

In order to make the amplitude A n (t) of the symmetric, odd modes corre- 
spond to the root-mean-square displacement over the length of the CNT (not 
in time) as: 



A n (t) = 




(1.4.17) 



the eigenfunctions t; n (x') have the following orthonormalization: 

/ U{z')^n{z')dz' = S mn . (1.4.18) 
Jo 

Inserting expression 1 1 . 4 . 1 7l into equation ll. 4.13[ multiplying by £, n (z'), integrat- 
ing over the length, and unsealing gives: 

Li-^ +i ^j An = u anj (1419) 

where a n = J Q t; n (z')dz'. The first four a n are: 0.83, 0, 0.36, 0. In the ideal case, 
even modes have a displacement that on average gives no change in capacitance 
and are therefore often not visible in experiments. 

Up to this point, a linear system was described. We now introduce nonlinear- 
ity, looking at two different contributions to the geometric nonlinearity of a CNT 
resonator. For the first contribution, the CNT is described as a beam [3"6l 157] , 
with no dc displacement or tension Ud c ,Td c = 0. For the second contribution, 
the two regimes from [38], i.e. the weak-bending regime and the strong-bending 
regime, are used to account for the dc displacement a CNT possesses and the 
tension it acquires from applying a gate voltage. 

Modeling the CNT as a beam, we see that when the amplitude of oscillation 
becomes larger, it is significantly longer at the extremal points of oscillation 
than at its equilibrium position. A segment that used to be of length dL = dz 
now has a length of dL = %/ dz 2 + du 2 w dz + \{jj^) 2 dz. The CNT becomes 
longer by an amount AL and this causes a change in tension of EAAL/L. The 
associated added tension T'* c is: 




(1.4.20) 
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where b n = J Q ^ ) 2 dz' , with the first four 6„ calculated to be 12.3, 46.1, 

98.9, and 171.6. The scaled restoring force per unit length due to the tension 
T'* c is given by: 

_ ft* _ rpl* v ^ac M A 2\\ 

which we add to equation ll.4.131 Multiplying the augmented equation by £ n (z'), 
integrating over the length and unsealing, yields the well-known Duffing equa- 
tion: 

^n^n ' abeam^n^-n — a ^nt {L.^L.ZZ) 



dt 2 Q dt ™ " -o^„^ n pA 



with 



b 2 E 

abeam = TTTi ■ (1.4.23) 



2L*pCn ' 

Here the identity J* ^g§r£,ndz = — J* {^f) 2 dz = —b n is used. The factor 

c n = f Q £ n (z') 4 dz' arises from our definition of a: if the analysis had been 
started not with the unsealed version of equation ll.4.21l but with f r L = —mabeamU 3 

the Duffing term would have become abeam ( Jo £n(z , ) *dz'j u 3 . The first four 
c n values are 1.85, 1.69, 1.64, and 1.61. From equation 11.4.231 it is clear that if 
the CNT resonator is modeled as a beam the nonlinearity is positive and thus 
always gives a stiffening spring behaviour. For our example CNT in table 11.11 
it results in a nonlinearity for the fundamental mode of 9.6 ■ 10 34 N/(kg m 3 ). 

Until now, we have assumed that the CNT acts as a beam with no static 
displacement and no static tension. In CNT devices, these assumptions are 
often not valid: by applying a dc voltage to a nearby gate, the CNT can ac- 
quire significant static displacements that can be larger than the CNT diameter. 
As described in [38], with no initial tension, there are two regimes describing 
the bending of the CNT: the weak-bending and the strong-bending regime. In 
the weak-bending regime, the CNT acts as a beam at low gate voltages, and 
both abeam a n d a wea k contribute to the geometric nonlinearity. It undergoes a 
transition to string-like behaviour at larger gate voltages as the CNT becomes 
strongly bent by the force from the gate. The crossover from the weak-bending 
to the strong-bending regime occurs when the restoring force due to tension is 
equal to the restoring force due to bending rigidity, Td c = EI/L 2 , or cquiva- 
lcntly, when the CNT displaces by a distance equal to its radius. At this point 
the contribution abeam to the geometric nonlinearity, arising from the tension 
induced by the oscillation of the CNT, can no longer be described by equation 
11.4.231 This is because, at non-zero dc deflection in the strong-bending regime, 
the oscillation of the CNT does not lead to increased tension when moving 
away from the gate. The gate voltage V* at which the crossover from the weak- 
bending regime to the strong-bending regime occurs is found from l' dc w 871 
and for the example CNT in table 1 has a value of V* = 2.1V. 

At low gate voltages, the CNT is in the weak-bending regime with Td c -C 
EI/L 2 . In this range of gate voltages, the restoring force is dominated by the 
bending rigidity, leading to mechanical modes with resonance frequencies at 
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ratios close to that of an ideal beam. Starting from a straight CNT (i.e. no dc 
deflection at zero gate voltage), the tension induced by the gate is [38] : 

T ^ ak = 6M80EP F ^ (1A24) 
and the mechanical resonance frequency of the lowest mode is given by: 

[EI Pl T 

'^M + 07 ^ f^V (1 4 25) 

For a CNT that begins with no deflection, the resonance frequency at low gate 
voltage in this beam-like, weak-bending regime scales with u> oc Vg. Note that 
this quartic gate voltage dependence arises due to the fact that the electrostatic 
force from the gate is perpendicular to the axis of the CNT, and therefore does 
not begin to add tension until the CNT develops a static displacement. If the 
CNT already has a static displacement at zero gate voltage, due, for example, 
to buckling from compressive strain, or from a downwards curvature of the 
substrate at the edge of the trench, then the resonance frequency in the weak- 
bending regime shows a w oc V g 2 gate dependence. From expression 11.4.251 for 
the resonance frequency, wc obtain the gate voltage dependent spring constant: 



EI$ in „ L 3 Aft fdCg 



We obtain the Duffing nonlinearity parameter by taking the second derivative 
of equation 11.4.261 with respect to displacement. We do this by using the fact 
that the displacement of the CNT resonator leads to a change in capacitance to 
the gate. This in turn is interpreted as an effective change in the gate voltage, 
since only the gate-induced charge q c = C g V g is relevant. 



1 d 2 k 1 (V a dC a \ 2 d 2 k 

a 



2m du 2 2m \C g du J dV 2 

(1-4.27) 



60480EI 2 C 2 p V du 
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For the example CNT in table 11.11 this results in a nonlinearity of 8.7 ■ 10 
(Vg in Volt) 4 N/(kg m 3 ). Note that this gate-induced contribution to the non- 
linearity in the weak- bending regime is much smaller than the contribution from 
the intrinsic gate- voltage independent beam nonlinearity (equation II .4. 23| ). and 
thus the nonlinearity in the weak-bending regime is dominated by that described 
in eauation ll.4.231 

As the gate voltage is increased further, the CNT enters the strong-bending 
regime, in which its static displacement exceeds its diameter, and the restoring 
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force from the induced tension exceeds the bending rigidity. In this regime, the 
mechanical modes begin to behave as those of a string under tension, approach- 
ing a harmonic spectrum at large gate voltages. In this regime, the gate induced 
tension is given by |38j : 

t EA 

with a resonance frequency for the lowest mode given by: 

L 

P A -L^A\\96 du/^- (1A29) 

2 

In this strong-bending regime, the resonance frequency scales as w oc V g 3 . The 
gate voltage dependent spring constant is given by: 

L \ 96 J \ du 
and we obtain the following for the Duffing nonlinearity parameter 




l T dT° ng _ 71 ( [EAdCg 



k strong = T {-) ^) Vf, (1.4.30) 



^■strong 



SmLCl {du ) [96j Vg ■ (LA - iL) 

For the example CNT in table ITTTl this results in a nonlinearity of 1.9 • 10 28 • 
(V g in Volt) I N/(kg m 3 ). 

Until now, we have assumed that the CNT includes no initial tension, i.e. 
that its relaxed length is the same as the length of the trench. For positive initial 
tension, corresponding to a tensile stress, the CNT behaves already as a string at 
zero gate voltage, and shows a fiat resonance frequency vs. gate response. This 
flat gate response of the resonance frequency continues until the gate induced 
tension exceeds this initial tension and begins to tune the resonance frequency, 
as shown in figure ll~4~3"l 

If the CNT relaxed length is longer than the distance across the trench, the 
CNT experiences a compressive force. For compressive forces below the Euler 
instability, this results in a reduction of the resonance frequency at low gate volt- 
ages. As the compressive force approaches the Euler instability, the resonance 
frequency at zero gate voltage vanishes, as shown in figure [l~4.3l For compressive 
forces beyond the Euler instability, the CNT buckles and acquires a non-zero 
static displacement even at zero gate voltage. Beyond the buckling point, the 
resonance frequency becomes finite again, and the CNT exhibits bending-like 
modes around the buckled beam shape. In this case, the bending modes show 
an uj oc Vg gate dependence due to the non-zero initial displacement, as dis- 
cussed above. The modes of such buckled beams can display peculiar mechanical 
modes similar to those of a hanging swinging chain. A complete analysis of the 
mechanical modes of strongly buckled CNTs can be found in [39] . 
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In the strong bending regime, the gate voltage dependent spring constant 
given by equation 11.4.301 results in a significant quadratic contribution to the 
restoring force. This suggests that the quadratic renormalization of the cubic 
force coefficient may play an important role in this case, something that should 
be the subject of further investigation. 




2 4 6 

Gate voltage (V) 

Figure 1.4.3: Simulation showing the resonance frequency as a function of gate 
voltage for three values of residual tension Tq. At a high positive residual tension 
of Tq = 100 (dotted line), the resonance frequency only starts to increase when 
the tension due to the gate voltage is similar to the residual tension. At zero 
residual tension Tq — (solid line), the tension due to the gate voltage first has 
to overcome the bending rigidity in the weak-bending regime with luq oc Vg, 

2 

and then the strong-bending regime with luq oc V g is entered. At a negative 
residual tension Tq = —39.4 signifying residual compression (dashed line), the 
resonance frequency increases significantly with gate voltage even at low gate 
voltages. 



1.4.3 Electrostatic nonlinearity 

A second origine of nonlinearity comes from the electrostatic interaction between 
the CNT and the gate electrode. If for a moment we consider the CNT not to 
experience Coulomb blockade, but to be a perfect conductor, then it screens the 
voltage applied to the gate. The electrostatic force between the CNT and the 
gate was already given in equation II .4. 101 as: 



_ dU_ldC g 2 



In the case of a cylinder above a plate, the capacitance per unit length c g = C g /L 
is given by [38] : 

«*(*) = 2 ?° — : (1-4-33) 

arccoshl 



h g —u(z) 



23 



where h g is the distance between the CNT and the gate. For the example 
CNT in tableOthis gives values of 5.9 • 1CT 21 F/nm, 3.5 • 1CT 23 F/nm 2 , and 
2.9 • 10~ 27 F/nm for the first, second, and fourth derivative of the capacitance 
to gate with respect to displacement respectively. The force acting on the CNT 
due to electrostatics can be expanded as: 

F ES = ~Ak ES u = - (Ak ES . + ma ES u 2 + 0(u 4 )) u, (1.4.34) 



and adds to the spring constant as: 



dF ES 

du 
ld 2 C 



2 du 2 9 



V 2 (1.4.35) 



The electrostatic force thus gives a contribution to the spring constant that is 
opposite in sign to the contribution from the restoring force due to bending 
rigidity and tension. Note, in particular, that the electrostatic forces do not 
act like a restoring force: for example, if the CNT is displaced away from its 
equilibrium position towards the gate, the electrostatic force pulling it towards 
the gate increases. This has the effect of an "anti" -restoring force. In the 
case that this electrostatic "anti" -restoring force is larger than the restoring 
force from the mechanical deformation, the system can even become unstable, 
leading to a phenomenon of electrostatic "snap in" , in which the device is pulled 
uncontrollably towards the gate until it makes physical contact. 

In the case that the restoring force from mechanical deformation dominates, 
the system is stable, but the effective spring constant is reduced by the electro- 
static interaction. The reduction in the spring constant leads to a decrease in 
resonance frequency with gate voltage. When the CNT is close enough to the 
gate so that d 2 C g /du 2 is large, this decrease in resonance frequency is larger 
than the increase in resonance frequency with gate voltage due to increased 
tension. 

The nonlinearity parameter due to the electrostatic force is: 

1 d 3 F ES 
6m du 3 

= — ^fV 2 (1.4.36) 

12m du 4 9 y ' 

The electrostatic nonlinearity parameter is negative and always gives rise to a 
softening spring nonlinearity. For the example CNT in table[T7TJ the nonlinearity 
due to electrostatics is -3.2 • 10 28 • (V g in Volt) 2 N/(kg m 3 ). This is of the same 
order of magnitude as the geometric nonlinearities a wea k and ex strong-, but has 
an opposite sign. 



1.4.4 Nonlinearity due to single-electron tunneling 

In the previous section we described how single-electron tunneling causes a 
softening spring nonlinearity. Here, we describe how it is also an additional 
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source of nonlinearity. Writing the force on the CNT due to single electron 
tunneling as: 



Fset = -Ak SET u = ~ (Ak SET ,o + ma SET u 2 + 0(u 4 )) u, (1.4.37) 

where AksET is given by equation ll.3.5[ gives for the nonlinearity due to single- 
electron tunneling: 

a SET = ±fl^r = J_ <?Ak SET = J_ (V^dCA 2 <?Ak SET 
6m du 3 2m du 2 2m \C g du J dV g 

We now use the expression from eauation ll.3.5l and insert it into eauation ll.4.381 
Assuming that the gate voltage difference in which the charge transition occurs 
is much smaller than the gate voltage, and neglecting terms of d 2 C g /du 2 , we 
arrive at: 

asET = .JL ^-^r) (d^Vdm 
SET 2m Cf<7 tot V du J dV 3 ' 1 ^ j 

Equation 11.4.381 shows that the nonlinearity parameter is proportional to the 
curvature of the change in spring constant and consequently to the curvature 
of the change in resonance frequency. Figure 11.4.41 illustrates how the sign of 
<%set changes with gate voltage: in the middle of the frequency dip caused by 
single-electron tunneling the curvature of the change in resonance frequency is 
positive, leading to a positive nonlinearity, whereas on the left and right side of 
the frequency dip, both the curvature and the nonlinearity are negative. 

According to equation 11.4.391 the nonlinearity parameter depends on the 
third derivative of the average occupancy with respect to gate voltage. Figure 
ll.4.5l shows. using the values from the example CNT in tablc ll.ll how by chang- 
ing the gate voltage a few mV, the nonlinearity parameter changes dramatically 
and even changes sign. For Vcnt <C V g and a charge transition that occurs in 
AV^ = 10 mV the example CNT has a minimum (i.e. most negative) nonlinear- 
ity due to single-electron tunneling of a>s E T = —2.7 • 10 31 • (V g in Volt) 4 N/(kg 
m ). The maximum (i.e. most positive) nonlinearity is as E T = +8.0 ■ 10 31 • 
(V g in Volt) 4 N/(kg m 3 ). At the inflection point of the frequency dip, the cu- 
bic, Duffing nonlinearity due to single-electron tunneling vanishes. At this point 
other nonlinearitics, such as the quadratic nonlinearity, will dominate the nonlin- 
ear dynamics, which suggest further theoretical and experimental investigation 
of this regime. 

In table 11.21 the magnitude of the four different contributions to the non- 
linearity parameter for the example CNT in table 11.11 is depicted. The table 
demonstrates that the nonlinearity due to single-electron tunneling, a$ E T, is 
similar in size to the beam- like geometric nonlinearity, Ubeam- Which nonlin- 
earity dominates depends on the gate voltage and the gate voltage range over 
which the charge transition occurs. The nonlinearities due to weak-bending and 
strong-bending and the electrostatic nonlinearity have almost the same value, 
but opposite sign, and can be neglected for a quantum dot CNT resonator. 



25 




Figure 1.4.4: Comparison of two effects caused by single-electron tunneling: 
dips in the mechanical frequency, and mechanical nonlinearity. When the cur- 
vature of the change in resonance frequency is negative, so is the nonlinearity 
parameter. A positive curvature of the change in resonance frequency leads to 
a positive nonlinearity parameter. 

1.5 Parametric excitation and mode coupling in 
carbon nanotube resonators 

The previous section showed how geometry, electrostatics, and single-electron 
tunneling give rise to a Duffing-typc nonlinearity. In addition to a nonlinear 
response to large powers at the resonance frequency, nonlinear effects can also be 
observed when driving the resonator at frequencies different from its resonance 
frequency In the next subsection, we look at how nonlinearities in CNTs can 
result in parametric driving and amplification. In subsection ll.5.2l we consider 
other modes of oscillation, either in another dimension of oscillation, giving rise 
to nonplanar motion, or as higher order modes in the same direction. Through 
the nonlinear behaviour of the device, these other modes become coupled to the 
fundamental mode. As we will see, driving these other modes can generate, for 
example, extra tension in the device, which can then shift the frequency of the 
fundamental mode. 



1.5.1 Parametric excitation 

A resonator is parametrically excited when besides a driving force around ujq 
there is also a change in one of its parameters around 2loq. A well-known example 
entails a child on a swing which alters its center of mass during swinging. In 
this analysis we only consider a modulation of the spring constant, so that 
k = ko + k p cos(wpi). 
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Figure 1.4.5: Calculated average occupancy (left axis) and the nonlinearity 
parameter asET due to single-electron tunneling (right axis): olset depends on 
the third derivative of the average occupancy with respect to gate voltage, thus 
going from negative to positive to negative again in only a few mV. 



The pumping of energy into the motion through parametric excitation is 
maximum when done at twice the resonance frequency, but can be done at any 
frequency for which oj p = — uiq, becoming less effective for larger n. The equation 
of motion describing parametric excitation is: 

d 2 u u Q du 1 , . , F ac 

—r + — — H (fc + k p cos{uj p t)) u = cos(wi). (1.5.1) 

ot z Q at m m 

In the case of capacitivc driving at ui p , the change in the spring constant can be 
described as: 

Afc = -— = -^{(vt) 2 + 2V^cos{u J pt) + (V g ac fcos{2u J pt)) 

= Afc + kp cos(wpi) + k'p cos(2w p £), (1.5.2) 

yielding an extra modulation of the spring constant k' p at 2tu p . For low drive 
powers, Vg C <C V ff dc , so that k' p k p , and k' p can be neglected. The most 
effective parametric excitation can then be achieved by driving at lo p = 2u]q. 

Recently, parametric excitation has been applied to a CNT resonator [40] . A 
CNT resonator is especially suited for parametric excitation because the change 
in spring constant with drive voltage, dk/dV g , is large (see equations II .4.261 and 

1.5.2 Mode coupling 

Until now we have only considered oscillations of the fundamental mode taking 
place in the x-direction. In this section we discuss coupling between the fun- 
damental mode in the x-direction to the fundamental mode in the y-direction 
[411 1421 112] and between the fundamental mode in the x-direction and higher 
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Table 1.2: Comparison of the magnitude of the parameters for the four contri- 
butions to the nonlinearity parameter for the example CNT in table ll.ll For the 
nonlincaritics due to single-electron tunneling the gate voltage range AVg over 
which the charge transition occurs is indicated. The nonlinearities in the weak- 
bending, beam-like (strong-bending) regime are not calculated above (below) 
the crossover voltage V* . 



modes in the x-direction [43] . Both types of mode coupling can generically be 
described by the following two coupled equations: 

, d 2 u du „9 4 u . . d 2 u 

pA W + v m + D d^ ~ {Tdc + Tu + Tv) d^ = fu > (L5 - 3) 

, d 2 v dv d A v , , d 2 v 

pA W + + D o^ ~ (Tdc + + T ^ = (L5 ' 4) 

Analogously to equation II. 4.141 the displacement of the other mode is denoted 
by: 

v ac = B{t)i>(z). (1.5.5) 

Here, B(t) is the slowly varying amplitude of the other mode, root-mean-squared 
along the CNT, and tp( z ) depicts the modeshapc. 

The different parts of the tension on the CNT by the oscillations u and v 
are given by: 

t - = §z£(^ mM ) 2 ^ <l5 - 6 > 

T ° - iff (^""H"*- (L5 - 7) 

Here T u causes the nonlinearity for a beam as described in section ITT41 
The tension arising from the oscillation of v is: 

r . . /:{%)'<■ 

= ^B* + ^&cos(2u> v t), (1.5.8) 
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remembering J Q L ^^7^ dz — b^/ L. This means there is both a contribution to 
the dc tension and there is change in the ac tension at 2ui v . 

Let us first look at the contribution to the dc tension in the weak-bending 
limit. Using equation 11.4.251 and 11.5.81 in the limit of T v <C EI/L 2 the change 
in resonance frequency to mode coupling is: 



Au;^ k = J — ^B 2 . (1.5.9) 

For CNTs, the high Young's modulus, low mass density, high ratio of A/ 1 = 
4/r 2 , and small length cause a large change in resonance frequency due to 
nonplanar motion. For the example CNT from table 11.11 this gives a tuning of 
16 MHz/nm 2 when the other mode is the fundamental mode in the y-direction. 
For the third mode the tuning of the fundamental mode for the example CNT 
is 127 MHz/nm 2 , arising from a higher V n . 

In the strong-bending limit with T v <C J 1 **™™ 3 : using equation 11.4.291 leads 
to the change in resonance frequency due to mode coupling: 



Aw 



strong 



Tr 



on a 

* dc 



4 • 2^L3 I du 



V g 3 B 2 . (1.5.10) 



As the tension due to the gate electrode increases the added tension from mode 
coupling is less important. At V g — 5V the tuning of the fundamental mode of 
the example CNT due to another fundamental mode and due to a third mode 
is reduced to 3.3 MHz/nm 2 and 26 MHz/nm 2 respectively. 



1.6 Conclusions 

Carbon nanotubes display a rich variety in nonlinear effects. Because of the 
importance of the measurement apparatus on nonlinear effects, we started this 
chapter by displaying the multitude of methods to detect the flexural motion 
of a CNT resonator. In the next section, we illustrated how single-electron 
tunneling in quantum dot CNT resonators, combined with their low mass, leads 
to pronounced softening-spring behaviour and significant damping. In section 
four, we covered the four different contributions to the nonlinear oscillation of a 
CNT resonator: the beam-like nonlinearity, the nonlinearity due to gate-induced 
tension, the electrostatic nonlinearity, and the nonlinearity due to single-electron 
tunneling, the combination of which is unique to CNT resonators. In the final 
section we show how the large response of the resonance frequency of a CNT 
resonator to a change in gate voltage or tension makes CNT resonators ideally 
suited for parametric excitation and mode coupling. 
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